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Abstract 


Direct electrical conductivity of the ionosphere is examined 
mainly concerning the effect of Hall current. It is suggested that the 
effect of the electric field produced by vertical component of the Hail 
current, together with the smaller effect by its horizontal component, 
may construct, along the magnetic equator in the ionosphere, a narrow 
region, in which the conductivity near the level of the E layer is much 
greater than in the other latitudes. This result is considered to be 
very fit for the distribution with latitudes of the range of diurnal 
variation of the terrestrial magnetism. = 


1. General formulae. It is well known that the presence of a transverse magnetic 
field causes certain anomalies in the flow of electric currents in ionized gases. When 
an electromotive force is applied in the gases, the Hall current tends to flow, and 
that reduces the conductivity of the material. We shall here evaluate the magnitude : 
‘of the effects in the ionosphere. We represent charged particles and neutral particles, 
by ae 132: epootvely_2 and. let ¢1, ¢: be velocities of Bo Ni, Nx be oie 


SS the magnetic field respectively. 
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where p=pressure H=the earth’s magnetic field 
C=peculiar velocity, Ci=e:—c, C:=e:— 
o=e,H/m Ty =1/¥, =m me2nD12/kT 
n=number density =2+72 
Dy.=—diffusion coefficient for ion: we adopted the value by Sutherland’s model 
of inverse fifth law, considering the effect of electric charge [2]. In this 
calculation the molecular weights of neutral molecules and ions are both 
taken to be 20. 
é,=electric charge of particle 1 (emu): 
for positive ion e,=e, for electron e;=—eé 
and de £ ( a )+ “ede 2) 2 a red Fi Es) peso x A) (3) 
‘in which mF, mF, are external forces acting on charged and neutral 
rarticles, and dj.=di+d"12, dix and dz are components of diz perpendicular 
and paralle] to the magnetic field respectively. 
In the £ and Fl region, according to numerical calculation, the first and the second 
term of (3) is small compared with the other terms, and in the third term, as the 
effect of the gravity does not appear, we may take only the effect of the electric field, 


therefore dy, = — 2 (EO + 60x A) pe ele ; (4) 


where EH represents electric field. Let H=E®+ce.xH- 
then. (3) becomes _ . 
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Consideration on the Hail current. a 3 pa = = os = sigices 2 
‘From (6), charged particles ee to move not ae o Ae direction of EB, he 
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H.C. cannot flow, counteracted by resulting electric field, the conductivity takes the 
same value as if the transverse magnetic field is absent. But in the case of the 
ionosphere, as we cannot consider an ideal case that charged particles do not flow at 
all in the direction of HH, we shall treat the problem more precisely. 


3. Equation of conductivity. 

Let right-handed rectangular axis o-xys be taken, so placed that ox, oy, oz are 
directed to the south, east and upward respectively and the unit vectors in the direc- 
tion of ox, oy, oz and A be i, j, k and HM) respectively and let R=jxM, dip of 
geomagnetic field be ¢. 

Let number densities of electrons, ions, +-ions, —ions, be 2, 7%; N+, N— respec- 
tively and let n_=An., nmy=(A+1) ne, mi =(2A4+1)n. = Kn. (6) 

In the following calculation it is assumed that the axis of rotation and the 
magnetism of the earth coincide, the following two cases (a) and (b) are considered. 
(a) Only vertical component of the H.C. does not flow. 

Here we discuss conductivity of the ionosphere in middle and low latitudes 
only. © ; : 

(1) Conductivity of the ionosphere in the direction of E-W. 

When an electromotive force (O, Ey,,O), parallel to the magnetic equator, is 
imposed, then vertical component of the H.C. produces an electric Polereuoe which 
produces an electric field E.. 

Bc E= s qe 3) Substituting this relation in (5) we Ret 
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= | From ¢ these relations and (5) we have. B= H iy Es ar aby ee = 
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pi=(1/m){vi/(vE+o2)} etc, pi=pert KPa (12) 
on the other hand the condnctivity by ions only is, 
Oiy~e2(KNe/ i) {¥i/ (Ve +.02)} (13) 
(2) Conductivity of the ionosphere in the direction of N-S 
Under the assumption that there is no vertical electric current, we have 
Ox=Ne? € © Pipo/(pi © Cos? +P sin? ¢) (14) 
where = Po=1/m, po= Poot kKpin | (15) 


When an electro-motive force #=(E,, E,, O) is imposed, the electric current 

I= (1, Iy, O) is as follows 

T,=62° Exton? Ey, — Ty=dyx° Ex+oy* Ey (16) 

where Ony = —Oyx= —Nel*qu ° sin d Ss as 

is obtained in the same way. 

(b) The H.C. does not flow in any direction. x 4 
(1) Conductivity of E-W , a 
Let an electro-motive force be Ho=(O, Eo,, O) and E’ =(E, Oz E) be an electric 

field by polarization of the H.C., then H=(E&;z, Evy, Ez) 

Substituting this relation in (5), we get, for sCHE. 


G=% on LEO ae 2 4(- E, sin ¢+E: cos $c) R i By—on(—Ees sin +E. etd wed 


: —t(Ee cos $+, sin 9) 

then Wis od ssEu-t bibs ay UsediEct SE.+uE,) ete. 8) 

As, here, both horizontal and eric ERE of the HC; $32. 3 
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where =the number of electrons primarily produced in cm'/sec 
a=the recombination coefficient for electrons in cni’/sec 
8=the attachment coefficient for ions in cm3/sec 
o=the coefficient of photo detachment in 1/sec 
y=the recombination coefficient for ions in cm3/sec 
m=the number density of atomic oxygen 
W,=the number density of particles, for which electrons can attach, that is 
Na=No OF Ne=No2 
From above equations, we have n_=An, 
where A=(ané+— a oe +—;+ 2 + Bnctta—1)/(Oneno + 0%e— ane) (21) 
According to Penndorf [4], the transition layer is situated near the height of n=5el0!3, 
Oxygen is almost atomic in higher region than that, and almost molecular in the lower. 
According to Bates and Massey [5] in the F2 and E layer, 4=8-10-, 9.10 
respectively in the day time. 
Below the transition layer, neglecting dT/dt, di/dt, from (21) we have > 
A= Bn. © 2a/(0Ne+7NE +9) . | (22) 
where q=1/(1+2) | 
- p=0.14 (almost independent of the seattle angle of the sun) 
7=107~10° 
and B=4.3x 10-8 e107 where z ‘is height from 60km level measured in 
Swit or au ee 
As q is not greater than see we can 
neglect q and therefore ——_ 
ds EIXIO XE HIOna/LA e) ae 
Ss ‘From the above we get Table 1. aye 
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When v<10" ; 2=103, «>10, the effect of Table 2 
H.C. dissappears. n Pr 
oes 102 1 
B 1033 il 
Case (3, a, 2) Cc 1013 20 
Table 3 —_— 
Gy /Me (emu) Ox [8% 
n | 103 10" n 103 10! 
x 2 x 2 
eee 1 10 ee 1 10° 
0° 2.16 «10-18 2.32 ¢ 10— 0° L.23'e108se 1.32 «10. 
102 9.7 «10-2 1.94.6 10-9 10°. 5.5 «abo 1.1 -10 


30° EO) elie tone oe (me 30° 29 ~ “4,83 


Cry/Me (indep of x) 


-—— 10" 1083 194 Page Ee? _ 
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0 0 i> (@ is the dip of 
10° 8.3 +102! 9.16 + 10-2", 8.756 10-21 ~ the geomagnetic field.) 
30° 2.38+10-2 2.636 10-20 2.51 + 10-9 sheds oN 
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Case (3,b,1) —%y /me is almost independent of $. 
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Table 6 


Using these results, we -calculated the conductivities 


Height (km.) n of the E and F2 layers which have Chapman distribu- 
110 2.41 +102 tion, at midday at the equator. In the E and F layer, 
100 8.05 e 102 we assumed that at the level of maximum ionization, 
90 6.98 « 1038 n=6e 10", 4.310" respectively and the temperature 
80 4.82 « 10" is 300°K, 1,000°K ; local scale height is 10° cm, 5.42 « 10° 
70 2.05 « 10% 


cm respectively. Then the conductivity integrated 
with height is 1.6 10-2.) for the E layer and 1.33 10-"n,) for the F2 layer, where 
Neo is Maximum electron density for each layer. 


6. Examination of the result. 

We examine the electrical conductivity of the ionosphere mainly concerning the 
diurnal variation of the terrestrial magnestism. At first we consider that electron 
density is everywhere constant at the same altitude. 

As it is most probable that Sq-current flows mainly in the E layer or in and 

blow the £ layer, the following two cases are considered. 

(a) Only vertical component of the H.C. does not flow. 

If «<50 at the height of 2=10”, 10'3, and if e<100 at the height of n=10™, dy is 

considerably great near the magnetic equator and rapidly decreases at 10° of latitudes 

towards higher latitudes. According to the examination of the previous section, 
_« seems to be in these ranges, if « is greater than these ranges, this tendency dis- 

appears. At low latitudes o,, is small compared. with o,. Therefore, when an electro- 

motive force, which is nearly parallel to the equator, as in ze field near mutoaty is 

imposed there, J, is approximately determined by A: e Evy. a : 


ay bib’ a bid Adee eile 


(b) The H.C. does not flow in any direction. 
Below the height of ~=10¥, ao is es greater than o;, and.o, is almost 


_ independent of latitude. _ = = 


_ The horizontal component of fe HC. Bee by the Sa current i in n the day time: 


According to the statistical research by Prof. 

IN M. Hasegawa and Dr. M. Ota, as shown in Fig. 3, 
there is a narrow region along the geomagnetic 

B equator, about 20° of latitude wide, in which the range 

\ of the diurnal variation of terrestrial magnetism is 

SUN much greater than in the other region. On the other 

Fig.4 Idealized. Sq Garent System hand radio exploration of the E layer seems to show 
that the electron density of the Jayer approximately 

follow Chapman distribution, therefore there is no particularly great electron density 
in the narrow region along the geomagnetic equator. When we try to explain this 
distribution of Sq by the dynamo theory, it seems very reasonable to consider that 
the mode of electric currents in the E layer or in its vicinity i. in the region of 
10%<n<10", is between (a) and (b). sits 
In the real tenosphere, near the level of the E layer, there will be a narrow 
region, along the magnetic equator at that level, in which the putes ses is manic 
- greator than in the other latitudes. ~ : 


Moreover, it is considered appropriate that AS least in the daytime Kis of order 


= near the level = n= 10", ate 
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Circulatory Motions in the Ionospheric 
Atmosphere and their Relation to the S 
Field of the Terrestrial Magnetism. III 


By Sadami MATSUSHITA 


Geophysical Institute, Kyoto University 


Abstract 


The relations between the shifting of ionic clouds of the sporadic- 
E and the type of the daily variation of the terrestrial magnetism are 
discussed in more detail than in the first report. Moreover, a research 
is being made regarding the shifting of ionic clouds of the F2 layer. 
From these results, it becomes necessary to consider the semi-. 
diurnal variation of the ionosphere. Therefore, the semi-diurnal 
variation and its relation to the S field of the terrestrial magnetism 
are discussed. Also a research of the dynamo effect due to both the 
¥: ctigey and the semi-diurnal variation is being conducted. 
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Daily Variation of the Terrestrial Magnetism . : 


—— 


‘The relations between ‘the shifting of ionic clouds of the sporadie-E &s) and | 
t type af the daily, variation, of is terrestrial magnetism have been discussed _ 
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Fig. 1 Classification of the type of the daily variation 
of the terrestrial magnetism. p, M and E are 


2 
? 2 owing to Hasegawa’s classification (2). Other 
= notations are as follows: 
a p! when maximum is in the forenoon. 
g pu when maximum is in the afternoon. 
p! s p (M) when difference from maximum to mini- 
3 mum is small. 
F M (p) when difference from maximum to mini- 
mum is smaller than p(M). 
g M(E) when variation is small in the forenoon 
? ww but increases in the afternoon. 
= Types with suffix s indicate some disturbed cases 
with somewhat intense decrease at night. 
3 
PM) = 
= 
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In Fig. 2, the relations between the shifting and the type, mentioned in the first 
report, are also obvious in this period. Namely, in the case of the P type and 
especially the Ps type, the change of direction of the shifting for the four time-groups 
is conspicuous, though it is obscure in the case of the M type and the storm time. 

From the aspect of these changes of the direction of the shifting for the case of 
the P type, it becomes.necessary to consider the semi-diurnal variation of the ionosphere. 
Tt will be discussed in paragraph three. 

In addition, frequencies and the time for occurrence of the intense and sharp 
decrease of the value of the horizontal component of the terrestrial magnetism in the 
case of the Ps type adopted to this statistical 
study during June solstice and equinoxes of 1949 
are shown in Fig. 3. It occurs most frequently 
at 19h., and this may have some relation to the 
Sp_ field. The obvious southward shifting at 
night of the Ps type may be an evidence of the 
effect of the Sp field, beside the effect of the 


; 

; semi-diurnal variation. Fig. 4 is frequencies of 
: TT Figen 0 BRS ee 0.510715 20 
: JST each type of the 


Fig. 3 daily variation 
of the terrestrial magnetism, Adopted to this statistical 
study during June solstice of 1949. 

The simultaneous variation of the Es region 

ionization in all observatories of Japan was-discussed — 
in the first report. Frequencies and the time of occur- 
3 rence of the simultaneous variation during June solstice 
of 1949 are obtained as in Fig. 5. It is interesting to_ 
note from the os ae the Seas variation = 


Frequencies of’ occurrence of the 
simultaneous variation of the Es in 
the four time-groups for each type of 
the daily variation of the terrestrial 
magnetism are shown in Fig. 6. It 1s 
interesting, -in the figure, that the 
simultaneous variation seems to occur 
frequently at the M type, and that 

[TEV [LUM ID0N IT0N the variation does not occur from the 
Fig. 6 Frequencies (in percentage) of occur- sun-set to the sun-rise, except the 


rence of the simultaneous variation of case of the storm time of the ter- E 
the Es region ionization in the four : 
time-groups for each case of the P, the 
Ps and the M type of the daily vari- 
ation and the storm time of the ter- 
restrial magnetism during June solstice 
of 1949. Dotted line shows the cases 
when the variation is somewhat in- 
conspicuous. When the variation of the criti- 

_ 1: 01-07 h. II: 07-13 h. cal frequency of the F2 layer (/°F2) is 

Til: 13-19 h. IV: 19-01 h. 


restrial magnetism. 
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2. Shifting of Ionic Clouds of © 
the F2 Layer 


investigated, some increment of f°F at 


midnight, in spite of no change 
of the height, is noted. The q 
time of such occurrence is from sy 
19 to 03h. Burkard (3) called _ 
‘it ‘ sporadic - F layer’. When 4 


| such variations at the five iono- 
_ spheric observatories in J, ust: ; 


JUNE 4, 1949 
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or the north, as in the case of ionic clouds of the Es region ionization. The time and 


frequencies of occurrence of such shifting in June solstice and equinoxes of 1949 are 
shown in Fig. 8. 
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: 19 21 23 1 3 JST. 19 aes) 3 
F Fig. 8 Times (abscissa) and frequencies (ordinate) of occurrence of the southward 
; shifting (the upper part) and the northward shifting (the lower part) of ionic 
i clouds of the F2 layer at four observatories (Yamagawa is omitted, by the 


reason of few occurrence) in June solstice (left) and equinoxes (right) of 1949. 
The unit scale of ordinate is 2. Black squares are conspicuous cases of the 
- shifting, and white squares are _ the cases that are somewhat inconspicuous, 

W: Wakkanai, se Fukaura, S: Shibata, ‘K: Kokubunji._ 


_ The apparent Jase of the shifting is — 1080 to 360. em/h, anid te 540 jm/h 

as the mean of all. (cf. 370 km/h as the mean, in the case of the Es.) Most recently, 
the speed of sporadic clouds was obtained from the observation at Baker Lake by 
Meek (4). It was from 600 to 1300 km/h at the height from 250. to. 300 km. He 
adopted 1200 km/h as the mean speed at the height of 300 km. It nearly corresponds 

> the abovementioned result, coe 400 km/h obtained by him for ga 
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Numbers with question marks show frequencies for the cases when the shifting is not 
so obvious. As the examples are quite few, it is difficult to obtain accurate conclusions 
from this result. However, the southward shifting tends to occur at night for the P 
and Ps type. It resembles the case of the Es. F urther, during storm time, the shifting 
—both the southward and the northward—occurs somewhat clearly, and in the case 


of the M type it does not occur 


F2 Layer at- ale 
Es Region S. Shifting N. Shifting When the southward or 
June ~ lpouinoxes June . [Equinoxes northward shifting of ionic clouds 
Solstice Solstice 

of the F2 layer occurs, frequen- 

S. Shifting 4(3?) 3 1(?) 0 - ; 
N. Shifting 1 0 1 9 ° cies of correspondence between 
Simultaneous | 3(2?) 0 0 0 the shifting of the F2 layer and 
Not Occur 2(1?) 11(5?) 4 1 that of the Es region in the same 


night are shown in the left table. 


Though the cases are quite few in this table, it seems the southward shifting in the 
F2 layer corresponds the sameward shifting in the Es region, when the shifting occurs 


in both regions. 


3. On the Semi-Diurnal Variation — BSS 3 
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dine the second report 0, Ts diurnal variation of ce ionosphere due to circu-. 
latory motions was discussed and the dynamo action owing ‘to the variation ‘was 
calculated. But, from the aspect of the abovementioned change of the direction of 
the shifting of ionic clouds of the Es for the case of the P type of the daily variation | 
of the. terrestrial ‘magnetism, its becomes: necessary to consider _ the 


= variation of the Enos Se ee tee ee 
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Fig. 9 upper part: Méeridional distribution of the semi-amplitude of the semi-diurnal 
variation of the virtual height of the F2 layer in Mar. (left); June (middle) 
and Sep. (right) 1947. The horizontal lines indicate C, and the curves are 
C+DP%, (The curve in June is treated only in the northern hemisphere and _ it 
is transfered 15° to the north from the form of C+DPj.) Dotted curves are 


A+BP: of the diurnal variation (see-the second report). 
lower part: The phase in hours, at which maximum of the semi-diurnal vari- 
ation of the virtual height of the F2 layer occurs, in Mar. (left), June (middle) 


and Sep. (right) 1947. 
Abscissa: geographical latitudes 


And the amplitude of the semi-diurnal variation for equinox is nearly of the same 
order as compared to that of the diurnal variation. But the amplitude of the diurnal 
variation for June in the northern hemisphere is about three times larger than that 
of the semi-diurnal variation, as the mean. Accordingly, in equinoxes and in middle 

latitudes of the northern hemisphere in June solstice, the semi-diurnal variation must 
be given important consideration. This is favourable to the foregoing results of the 


shifting of the ionic clouds. 

In middle latitudes, the phase in hours is nearly same (about 13 h.) both in 

equinox and in solstice (see Fig. 9), though it has great difference in the case of the 
iurnal variation. From these results and from the point of view of the cause of the 

semi-diurnal variation in the upper atmosphere (5) (6), this variation may be caused 

‘by a tidal motion. But the form of this tidal motion will not be such a simple form 

as P:, but may be the form of Ps or the combined form of Pg and Pi. 


4 If it will be possible to assume that the phase for the semi-diurnal tidal variation 
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(for example, the D layer or the lower 


of the height of a layer lower than the E layer 
as compared to that of the E 


part of the E layer) has such a phase difference of 180°, 
as in the assumption for the circulatory motions, and that the form 


and the F layer, 
r is represented by the form of 


of the semi-diurnal tidal variation in the lower laye 
P2?—3P%, then it will be calculated by the dynamo theory (7), that the lower layer may 
produce Sq in spite of small electric conductivity, when the wind velocity of the layer 
is of order 100 m/s in middle latitudes, as in the case of the shifting of ionic clouds 
of the sporadic E region, and when the variation of the height is assumed to be the 


Ese EN ee Te en 


vertical rise and fall of an isobaric surface. 
If the relative pressure oscillation for the semi- -diurnal tidal circulation in the 
lower layer is assumed to be 
0.15 (P§—3P%) sin (2t+60°), 
there is a source of the air flow in the middle latitudes of the lower layer at 10 h. 
of the local time. The phase angle 60°—from this the abovementioned 10 h. is obtained 
—is determined by the assumption that the semi-diurnal variation in the lower layer 
has a minimum at 13 h., from the results that the semi-diurnal variation of the height _¥ 
of the F2 layer in middle latitudes has a maximum at about 13 h. as in Fig. 9._ 
From the assumption, the velocity potential is 
| 1.28 10 (P§—3P4) sin (2t—30°). _- 
ae Fig. 10 indicates the air-velocities for this velocity potential. In the EK. and -F layers, 
—— the direction of the wind is reverse to it in this figure. _ Accordingly, the direction of — 
eo the wind in the E and F layers at the north of the source or the sink of the air-flow 
= corresponds to the statistical results of the shifting of ionic clouds of the Es. = 
When the electric conductivity K is constant, the current t function may be | 
Phe (0. 38 P?-0.61 Peas P§) x 108K sin (2t—30° Sse es on ates 
= If the total electric current for this current function 3 is of order 62,000 ai 
= —— ne of ss is ee a : 
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In this case, the electric conductivity is calculated with the conception for only 
one current-seat, and so it will be necessary to consider the effect of the reverse 
electric current in the upper layer. If the electric conductivity of the layer for Sq 
is of order 10—-° e.m.u., assuming the wind velocity in middle latitudes of the layer is 
about 100 m/s (when the electric conductivity of the layer is much more, it will be 
sufficient with smaller wind velocity) and assuming the electric conductivity of the F2 
layer is of order 10-° em.u. (6), then the lower layer will produce the Sq field 


sufficiently, under the conception of the effect of two current-seats. 


4. The S Field due to both the Diurnal and the 


Semi-Diurnal Variation of the Ionosphere 


The diurnal variation of the ionosphere due to circulatory motions and the 
semi-diurnal variation owing to the tidal circulation have been considered separately, 
and the dynamo action have been calculated respectively. Now the dynamo action due 
to both the diurnal and the semi-diurnal variation will become necessary to consider. 

From the comparison of the amplitude of the diurnal and the semi-diurnal 
variation of the ionosphere, it may be estimated as the mean state of a year that the 
dynamo effect due to the semi-diurnal variation is 0.6 times larger than that due to 
the diurnal variation. Because, from the abovementioned results that the amplitude 
of the diurnal variation of the F2 layer is three times larger* than that of the semi- 
diurnal variation in June solstice, but both are nearly equal in Dec. solstice and 
equinoxes, the ratio of the amplitude of the diurnal variation to that of the semi-diurnal 
variation as the mean state of a year is 5: 3, ie. 1:0.6. Accordingly, the Sq field due 
to the two variations is calculable by the assumption that the foregoing dynamo effect 
for each variation contributes respectively to the production of Sq with such a dynamo 
effect as this ratio. 

Fig. 11 is the electric currents obtained from these estimations, when the electric 


conductivity at this time is about one-half of the conductivity in the case of the 


calculation of the electric curre 
to the electric currents obtained from geomagnetical data (7) (8). 
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nts for either one variation. It resembles generally 
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i= * In Japan, the amplitude of the diurnal variation in June solstice is not always larger 
It corresponds with the semi-diurnal change of the 
But in geomagnetically high latitudes, 
of the 


M4 than that of the semi-diurnal variation. 
direction of the shifting of ionic clouds of the Es. 
of the diurnal variation in June solstice is very larger than that 


| 
= 
5 the amplitude 
- gemi-diurnal variation. This will be discussed in near future. 
: 
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Fig. 11 The electric currents—unit, 1,000 amp.—associated with both the diurnal cir- 


culatory motions and the semi-diurnal tidal circulation. (shown for the northern © 
i soe 


epee ise the case of the Sp field too, both ches irae and the semi-iurnal variation: 
must be considered. The principle for the calculation is already mentioned in the 
= Ese 
= second report. — See : Se ee Ay eee an ae DP Ae se a 
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were analyzed in each season, and the meridional distributions of their amplitudes 
and phases were studied, using the data of the virtual height from many ionospheric 
observatories over the world. From these results, it may be most possible to consider 
that the diurnal variation depends on circulatory motions owing to thermal expansion 
and that the semi-diurnal variation depends on tidal circulations. 


These assumptions may be convenient to explain a cause of the production of 


the S field of the terrestrial magnetism. Namely, if a layer lower than the E layer 


(for example, the D layer or the lower part of the E layer) has such wind that the 
direction is reverse to the wind supposed from both the shifting of ionic clouds in the 
Es region and the F2 layer and variations of the height of the ionosphere, then the 
lower layer may produce the Sq field and the E-F layer may produce the Sp field, in 


spite of their small electric conductivity (of order 10-8 e.m.u. in the lower layer, when 


the wind speed is of order 100 m/s), as the difference of the dynamo effect for two 
layers. . = 


ee 


Conclusively, a cause of the S field of the terrestrial magnetism was studied in 
some degree from the motions of the ionosphere, though these discussions were not 
always complete. 


In conclusion, she writer wishes to express his hearty thanks to Prof. M. 
Hasegawa, Dr. T. Nagata, Dr. M. Ota and Mr. Y. Nakata for their many valuable 
advices and encouragement throughout this study, and also to Miss M. Kato and Miss 
M. Kanayama for their assistance in numerical computations. SS 
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Magnetic Anomalies and the Corresponding 
Magnetic Centres. I 


By Tsuneji RIKITAKE 


Earthquake Research Institute 


Abstract 


The writer proposed a method for finding the position, intensity 4 
and direction of a underground magnetic dipole from the correspond- 
ing magnetic anomaly on the earth’s surface. In the first place, the 
distribution of the vertical force is expressed as follows ; 

Ze 25(_\gacien ty 
With the aid of some relations concerning spherical harmonic functions, 
then, the coefficients of the spherical harmonic expansion of the 
magnetic potential on a sphere that entirely contains the magnetic 
matter are determined. The radius of the sphere which is included 
in the expression as a parameter adjusted so as to make the higher 


a i 


degree terms minimum in the next place. 
Applying the method to the magnetic changes which aceermuaied 
the eruption of Miyakezima Island of 1940, a reasonable conclusion 


is obtained. 
* 


1. Seer odncion. Though it ae not posable to determine uniquely the subterranean 
mass distribution to which the magnetic anomalies on the earth’s: surface | are due, 
some aeioes have been a aa bys a number of investigators for finding e mass a 
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magnetism of igneous rocks,” the first degree terms in the spherical harmonic expansion 
of the magnetic potential are so taeee that the higher degree terms can be neglected 
in spite of complex arbitrary forms of rock specimens. Hence, it is easily presumed 
that the dipole terms are most important in the magnetic anomalies due to magnetic 
matter of which shapes do not deviate markedly from the sphere. Accordingly, the 
method will be applicable for any magnetic anomalies so far as the said condition is 
satisfied. Applying the method to actual anomalies, it will be able to determine the 
depth of the magnetic centres that roughly agree with the centres of the magnetic 


matter to which the anomalies are due. 


2. Theory. We consider a sphere which is sufficiently large so that it contains 
entirely the magnetic matter considered. Denoting the radius of the sphere by a, the 


magnetic potential is given by 
W=as By (4/1) Panda cosm¢+by sinmd¢) (r>a) ee 
n=1 m= 


where Pam denotes Neumann’s spherical surface harmonics. Taking a coordinate of 
which origin coincides with the centre of the sphere, where 6=0 axis agrees with z 
axis, and using the relations” : 


(—2)"n! 


“2r(n—m)! \ Se oie COs es sin wt cos mudu,) 


= 2 ys. 
1 Pym Sin mp 2S — cosatéysin -P! sin mudi) —— 


7" Pm COS MG = 


Re ea ee 


the typical term of (1) can be written as 


: ae Le eco yaa mec ees ae 
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Table I 
—e 0 1 2 Om 0 1 2 
if (5S Se a |S 
00 2n 0 0 E00 0 0 0 : 
™ 0 m 0 ai 0 0 0 | 
Was 0 0 0 Li 0 T 0 
fe x 0 > ah 0 0 0 
a 0 0 0 gz 0 0 es 
fn x 0 => Pe 0 0 0 
ae 0 an 0 £30 0 0 0 
oa 0 0 0 gy 0 =. 0 
a 0 a 0 ee 0 0 0 
fn 0 0 0 gin 0 — 0 
On the other hand, we can obtain er i from the observation or — observed 
distribution of the vertical intensity is expressed ae F 
: ee ae a, Vania biel oe 
In case of 7, k<3 and n<2, (4) gives the next spe SS equations for aie ihe 
ents of the splierical harmonic expansion. se = ay Be teeny Bee 
24;°+-3a=Z oo , : _ es a we 3 —— ; 
Bai! +4)/ 3 a2!= —adZ9/Ox, ee | iy 2 zs 4)» Socrerset ie “ied a 
Bhi F babe pias ee] ae 


ee ay 0+ as 5,/ 3.42 = == @O°Zyo/Ox?, | ts “ezine nena a 
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Ly = —@;!a29 + (ayao!+-ay!a.2 4 by'b.2))43 

L2= —by!a2" +-(ay°bo!— by ax? +- ay De") 4 3 

E£= (Loay+ Lia} + Lob, ')/4H)2 
Adjusting a or z so as to make z)=0, we can determine the parameter a or z. The 
depth of the magnetic centre is thus obtained. Further, with the aid of (7), the 
position, intensity, and the direction of the dipole can be calculated. 

In practice, the coefficients are determined by solving (6), while a remains as a 
parameter. The quantities included in (7) are calculated with these coefficients for 
various values of the parameter. The parameter must be determined by means of 
trial and error method so as to make 2o=0. 


3. Test of the method. In order to test the method, a simple case is studied as 
follows. In the case of the anomaly in the vertical intensity due toa magnetic dipole, 
which is burried at a depth of D and directs downwards, we have 


Pass = 2/ D3, OL /Ox? = OZ 9/ Oy? = 12/:D= 
O20 /Ox= OZ 0/Oy = OZ 9/Ox0y = OZ 9/023 = OZ 9/020 = OZ 0 /Ox0y? = OZ /Oy? = 0, 
where the moment of the dipole is taken to be unity. 
In that case, the solutions of (6) become 
odes: ESF 3a Se Se iecae 
a= 75 (5—“T : a Ca 


1 1 2 1 1 2 
a= = a=b;=h=b=0. 


Then it is obvious that the potential due to the second degree terms becomes minimum 
by taking a—D. Thus, as assumed before, the anomaly is interpreted as that due to 
a magnetic dipole situated at the depth of D. As a,°=1/D* and a'=b,'!=0, it is also 
proved that the moment must be equal to unity and the direction is downward. 


4. Actual example. As an actual example, the writer analyzed the magnetic changes 
that accompanied the eruption of Miyakezima Island in July, 1940. R. Takahasi and 
K. Hirano”: measured the changes in the vertical in- 
tensity during the period from July 20 to 29. From the 
repeated observations at 50 stations, they showed the 
distribution of the changes as reproduced in Fig. 1. 
The writer, in the first place, expressed the dis- 
tribution in the form 
Z=A+Bx+Cy+ Dx? + Exy+ Fy? + Gx + Hey + Ixy? + Jy, 
where the origin was suitably taken at the centre of the 
island. The values of Z at intersection points of series 
of parallel lines which were drawn on the figure at 
about 1.2 km interval were read off. - Using these values, 
the constants A, B, C,....were determined by means of 
least square. From these constants, Zo, 0Zo0/0x, OZoo/Oy, 


Opa eS 
kan 


The changes in the vertical 
force that accompanied the 


eruption of Volcano Miyake- ....were calculated. Then, solving (6), a,°, a}, bi, soteee 


jima.- Unit; y (After R. ae 
Takahasi and K. Hirano) were obtained. The value of @ that makes z=0 was 


2 


= 24 = 


calculated by means of trial and error method, being determined to be 2.9km. Mean- 
while, Takahasi and Hirano concluded that the changes were equivalent to those which 
might occur if a sphere, of which centre was situated at the depth of 3 km, lossed 
its magnetic susceptibility. Hence, the present result agreed well with their study. 
Further, it was found that the obtained- direction and intensity of the dipole also 


supported their conclusions. The position of the dipole thus obtained was marked on ~ 


the map with a small circle. 


5. Conclusions. The writer studied a method for finding the depth, direction, and 
intensity of an underground magnetic dipole from the anomaly in the vertical intensity 
observed at the earth’s surface. Taking into consideration that the dipole terms 
should be the most predominant in the anomaly due to magnetic matter, of which 
shape does not deviate markedly from the sphere, the method may be applicable for 
analyzing actual anomalies. The position of magnetic centre thus obtained will be 
roughly agree with the centre of magnetic matter which causes the anomaly. 


In conclusion, the writer wishes to express his cordial appreciations to Dr. T. 


— 


Nagata for his kind criticisms. 2 
- (Read: May 22, 1950) 
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Magnetic Anomalies and the Corresponding 
Magnetic Centres. II 


By Tsuneji RIKITAKE 


Earthquake Research Institute 


Abstract 


The method described in the first report is extended to cases in 
which magnetic anomalies are given in dip. 
aD changes in magnetic dip that accompanied the Shizuoka 
Earthquake of 1935 are studied as an example. 


—_ 


1. Introduction. In the first report? of this paper, the writer proposed a method by 
which the position of magnetic centres of magnetic anomalies were directly determined 


from the observation of the vertical intensity of magnetic field on the earth’s surface. 


The method will be extended here to cases in which magnetic anomalies are given in dip. 


2. Magnetic anomalies in X and Y. In a similar way with the treatment of . 
anomalies in Z, the partial derivatives of X and Y (the See a eastward com- 


- ponents) at x=y= 0 are calculated as follows ; 


(nt+j+1 
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Accordingly, we have the next relations for the cases 2>n 


—Y 3a! = Xoo ; 
—3ay9—6a29+ 1/3 a? =adXuo/Ox . 
VY 3 be? = a0Xoo/Oy ’ 
9a+15// 3 ae} =PO2X 0/02 , 
Shi) +5,/ 3 be! = OX y9/Ox0y ; a5 
3a, +5)/3 ae! = POX /Oy 
45,°a+135a:9—30//3 as? = aO°Xyo/Ox e 
—15)/ 3 be? = BOX /Ox0y , 
15a;9+ 45a." a B0*Xo9/Ox0y2 , } 
—15,/3 by? =BOXqo/Oy> | 
and —b'—y/ 3 be} as. Vag o | 
V 3 be? = a0 Yoo /Ox ; 
—3a;"—6a9—-\/ 3a? =a0Yu/dy , 
3b;'+5)/ 3 de! = 2°02 V0/Ox? 
3ay3+5)/3 ag} =@07Voo/Ox0y , 
9b,1+15)/ 3 be! =P0Y/dy  , oye 
—15)/ 3 be? =BRY/Ox* , 
15a,9+45a," - = BO? V/020y , 
15/3 be? = AO ¥o/Ox0y? , 


45a,°-+135a29+30,/3 3 a2= 2°03 Yoo /Oy3 e 

_ where it is assumed that a=z as before. 

By solving these simultaneous equations (3) and (4) Be eg NSS we can determine | 

a;°, a, ...., where it must be assumed that Xo, OX00/Ox, Fe Ge O¥o0/Ox, .... are 4 


already given from Taylor expansions of X and ee In order to determine the magne: 
tic centre, the parameter a should oe then so as ‘to ‘make the potentiz 


= 


ft (40)= 
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Combining (6) of the first paper with (3) of the present paper, we can easily 


get the next relation 


- (: 2a" + 3a") + 7 (ay'+} ay Qs!) = foo ’ \ 
i 5) l 1 5) ~ = 
Toray 3a4'+ 4% 3 az Aone (3ay9+ 6a —77 3 ar?) = a0 fyo/Ox ; 
Gb) fay’ 3 by) 3 bs vol 
a. Vo 02°)— x 3 bs? = a0 foo/Oy > 
1 Bie ee pa 1 — 
g(a, +30a2°—5,% 3 as") —(9ai" +15 “3 as") =POf)/O% , 
oe! 517 3 be? 1 JY Me ol by ee at a A , 
eee ere ee re 262) = @O"foo/Ox0y , 
s a (6) 
(12a, +304, +5; ra as!) (Ban +5y/3 az") =@O'fy/Oy?  , 
—-Z(H5a1+90,/3 a. 1)—— (4501? +135a; °— 30) 3 de®) = 0308 fy/Ox8 
1 >= Se 
yz (158i! +30 s be!) +15, f 3 by? = BO? foo /Ox20y , 
g(a +30,/3 as!) — = (15a," +-45a,) = GO foo /Ox0y?, 
1 = ce 
ira anes +904 8} by”) +15 a by? = 03 fo /Oy? 


The coefficients of the spherical hprmonic expansion are determined as the solutions 


of these simultaneous equations. 


4. Actual example. Ye 


Kato” observed repeatedly the changes in magnetic dip in 


the vicinity of the epicentral area after the Shizuoka Earthquake of July 11, 1935. 
The differences of the values obtained in the survey in Aug. from that in July are 


The changes in dip 
after the Shizuoka Earth- 
quake. (After Y. Kato) 


reproduced in Fig. 1. In like manner with the treatment 
in the case of Volcano Miyakesima in the first report, 46 
was read off from the map from which foo, Ofo/0x, 

are calculated by means of least square. F rom. the relation 
Z),=0, the parameter a was determined to be 4.2km. The 
position of the magnetic centre or dipole was marked with 
a small circle on the map. The north seeking end of the 
dipole directs toward E7°S. Its dip is also obtained to be 
7°. The moment amounts to 44x10" emu. The result 
agreed well with Kato’s study in which he determined the 
position, direction and intensity of the corresponding dipole 
by means of trial and error method. However, it is not 


clear why the changes in geomagnetism as obtained here took place in connexion with 


: 
; 


_ earthquake occurrence. 


oR eae 


ee ee 


6. Conclusion 
The writer showed in this paper that the position, direction and intensity of a 


subterranean magnetic dipole can be found directly from the anomalies in dip as well 
as in Z that was already studied in the first report. Applying the method to the 
anomalies in dip in the case of Shizuoka Earthquake of July, 1935, some interesting 


(Read: May 22, 1950) 
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results were obtained. 
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Magnetic Transition Points of Volcanic Rocks* 


By Takeshi NAGATA and Syun-iti AKIMOTO 


Geophysical Institute, Tokyo University 


The knowledge of the magnetic properties of rocks composing the earth’s crust 
is essentially significant in the studies of secular variation and local anomaly of 
geomagnetic field. Among various magnetic properties of rocks, the magnetic transi- 
tion points is the most fundamental one. Consequently the change in magnetic 
susceptibility with temperature of rocks has been examined by a few investigators. 
(1) (2) (3). In the present study, the thermal change of susceptibility of a large number 
of volcanic rocks was measured by means of a ballistic method in a weak magnetic 
field, several times as much as the geomagnetic field. The specimens examined were 
thirty in total and they are the ejecta from volcanoes Huzi, Amagi, Mihara, Usame, 
Hakone, and Taga, each of which was analysed chemically and petrologically by 
petrologists. These specimens were heated up to about 660°C and then cooled down 
to the room-temperature. During the heating and cooling processes, their susceptibility 
was measured every 20°C. 

The results of this experiments show that the susceptibility of all specimens 
becomes apparently zero at about 600°C, which coincides with the Curie-point of pure 
magnetite. The mode of change with temperature, however, is not simple. Its thermal 
change seems to be classified, generally speaking, into the four types such as given in 
Figs. 1—4, as the typical examples. In Fig. 1, the susceptibility in heating process 
increases grandually up to about 400°C and then decreases abruptly, and the change 
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No. 1. H=2.13 Oc. No. 2. _H=1.80 Oz, 
(No. 23) Huji (No. 55) Usami 
Fig. 1. Reversible ordinary type Fig. 2. Irreversible ordinary type 
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im cooling prowess follows almost the same curve: this change will be called “reversible 
ordinary type”. i Fig. 2, the change in heating process is similar to that in Fig. 1, 


Table 1 
REVESSISLE TYPE 


pene m 


while its cooling does not accord with the heating one: this change will be narneci 
“irreversible ordinary type”. On the other hand, in Fig. 3 and Fig, 4, the thermal 
change of susceptibility is stepwise, the former reversible and the latter irreversible: 
these types will be called respectively “reversible and irreversible extraordinary type”. 

As the magnetic transition points, two points on the temperatures aale were 
taken: the one is the temperature denoted by 9, at which ferromagnetion of rocks 


disappears, while the other is the temperature denoted by 9,, at which — op UT\= 


_ maximum. (4). These characteristic quantities of the measured D pectmens were 
given in Table 1, and Table 2, where z, Lmay, and T= maz deme respectively initial 
specific susceptibility at 0°C, maximum specific suxeptibility, and the temperature 
corresponding to the maximum specific susceptibility. Suffixes o and ¢ represent the 
ordinary and the extraordinary types. Reversible types were summarized in Table 1 
while irreversible ones in Table 2. i 
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Oosima ; 2.06 

038 Meizi-Taisyo lava Basalt | 1.06 | 
51 Usami Two-pyroxene-andesite 0.84 
56 Usami Augite-hypersthene-andesite 0.64 0.87 


Augite-olivine-basalt 


Taga 


IRREVERSIBLE TYPE 


Speci : : Heatin, 
tlie Locality Rock : e 

27 Sans Olivine-basalt 2.41104 

49 Usami Olivine-two-pyroxene-andesite 1.44 

50 Usami Aphyric-two-pyroxene- 2.06 
andesite 

52 Usami | Hypersthene-bearing 1.27 
olivine-augite-andesite ~ 

3 * Augite-bearing $ 

: SS olivine-hypersthene andesite of 

54 Usami Olivine-bearing = 1.43 
two-pyroxene andesite 

255 saci Olivine-bearing 


two-pyroxene andesite 


Hakone 
New somma 


| Hypersthene-augite-andesite 


Taga — 


Two-pyroxene olivine-basalt 


‘Tt seems clear that the maenene: es oo caret: upon ‘te chemical 
. constitution of the rock. In eee to see ae Bee goers — 6. of variors 
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210 590 420 550 310 2.69 10.57 


210 620 370 570 220 2.61 
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LETTER TO EDITOR 


An Estimation of Diurnal Variation of VHF Radio Waves. 


In this letter, the writer reports that the diurnal variation in the vertical 
distribution of the refractive index due to the variation in meteorogical conditions, is 
expected to be the main cause of the diurnal variation of VHF radio waves. 

In the lit-zone, the receiving electric field are mainly composed of the direct 
wave, the reflected wave from the ground and the ground wave, and are shown by 
the Bullington’s formula, 


f= 28in-4- +j[04+R)4+0—-R)A]e(® 2), 


where E, is the free-space field intensity. 

The variation in the vertical distribution of refractive index changes the wave 
path and reflection angle. For this reason, the receiving electric field changes. 

In actual cases, the vertical gradients of (#—1)x10°/km are often 50 in night 
and 40 in daytime respectiAely, in the lower atmosphere, where n is the refractive 
index. In this case, the ranges of diurnal variation are theoretically given as follows. 
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Range _ Polarization = Medium under path © 
5 ab 322 erica Land 
3. “Horizontal Land 
: 8 ee s Vertical ae Sea 
eer ee ee Horizontal Sea. 
a wave-length; = 5m | fore : a 
om eo distance; 110 km. i Se eee 


antenna height; 200m 


These range agrees with the experimental results in its” cuir: of magnitude, 
“In near future, the expe on 
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